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The Kondo effect discovered originally in connection 
with the resistivity minimum in dilute magnetic alloys 
turns out to be one of the most interesting many- 
body phenomena in condensed matter physics J2]. The 
localized spin-degenerate level embedded in the sea of 
conduction electrons may influence the properties of the 
system in a dramatic way. The key feature of the Kondo 
effect is the formation of the Kondo, or Abrikosov-Suhl, 
resonance near the Fermi level 0- This resonance is of 
crucial importance for so different fields as the physics 
of heavy- fermion systems [3114 and the electron trans- 
port through quantum dots jfj, |fj . Direct observations of 
the Kondo resonance at metal surfaces by the scanning 
tunneling microscopy method Q enhances the interest to 
the problem. 

In contrast to our detailed understanding of the equi- 
librium properties of the Kondo system, much less is 
known about it's real-time dynamics. The time devel- 
opment of the Kondo effect in quantum dots after the 
voltage change has been considered in Ref. H It was 
shown in Ref. 9; that the probe of the charge state of 
magnetic impurity leads to suppression of the Kondo res- 
onance. However, the dynamics of the latter process has 
not been investigated in detail. In this work we study the 
evolution of the Kondo system after the measurement of 
the impurity charge. 

We proceed with an exactly solvable large- Nf degener- 
ate Anderson model |lfl lllj in the infinite-?/ limit with 
the Hamiltonian 



H = P 



P, 



(1) 



where c kv , /„ are the Fermi operators for the conduc- 
tion band and /-electrons correspondingly (we consider 
the case where the /-level lies in the conduction band), 
V/ y/Nf is the hybridization parameter, i/ = 1,2,..., Nf 
is the "flavor" index and P is the projection operator into 
the space with n/ = Yl v ftfv < 2- We consider the case 



Nt — > oo. In this limit the ground state is [ToL ITU 

/ occ \ 



(2) 



where \Q > is the noninteracting vacuum, i.e. a filled 
Fermi sea of band electrons with Fermi energy e^. 

The wave-function after we probe charge localized at 
the /-level depends upon the results of the measurement. 
If we have found the hole the wave-function directly after 
the measurement is simply 



ip{0) = |0 > 



(3) 



Notice that in the accepted approximation the evolu- 
tion of the system after the measurement probing hole at 
the /-level coincides with the model of the initially non- 
interacting band electrons where we suddenly switched 
on hybridization with the /-level. 

In this limit we have, actually, a single-particle problem 
with the effective basis |0 >= |fi > and |jfc >= flc kv \^l > 
|lfj|. The matrix elements of the Hamiltonian in this 
basis are 



iT|0> = -e/|0>+ V k\ k > 

k<k F 

H\k> = - efc |fc>+y*|o>, 



(4) 



describing a discreet level with E — —tf embedded into 
continuum, presented by a band of a finite width with a 
bottom at E = —Ep. It is the finite density of states 
at the bottom of thus appearing band that is responsible 
for the formation of logarithmic divergences and, conse- 
quently, the Kondo energy scale in the problem. 
The wave-function can be presented as 



i>{t) = a{t)\f > + J2 b(k,t)\k>, 



(5) 



with the initial conditions a(0) = 1, b(k,0) = 0. For the 
amplitude to find electron at the /-level, straightforward 
algebra gives 



1 

a(t) = -7T-- 
2m 



(III! 



Ll> + €f — S(w) 



(6) 
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FIG. 1: Contour used to evaluate integral in Eq. ®. 



where 



eh = ]T 



k<k f 



\v k \ 2 ^ 

w + tk 



(7) 



The transition from the discrete level in the electron band 
is done by shifting integration contour by an infinitesimal 
value +is. Thus instead of TV poles at the real axis, 
corresponding to discrete electron levels we obtain two 
branch points: bottom of the band u = D and the Fermi 
energy u> = Ep = 0. 

In Eq.© the integration contour can be closed in the 
lower half-plane, so the integral is determined by the sin- 
gularities of the integrand in that half-plane. We can 
take the cuts along the line lo = —iy, oo > y > and 
u> = D — iy, co > y > 0. Thus 



a(t) = e-WRj + Icut, 

poles, j 



(8) 



where j labels the poles of the function (u> + tt — E(w)) 1 
given by the solutions of the equation 



= - e / + E 



k<Ckp 



+ £fe + is 



and the corresponding residues are 



kp 



(9) 



(10) 



where ak (u>) = V k * /(u> + e^). Notice, that for large t only 
the contribution of the real pole lo = —Ek survives, to 
give 



*(*) = 



E 



'Ek)\ 2 



JE K t 



(11) 



Now we present an alternative point of view at the 
process of tunneling into continuum, which will clarify 
the meaning of Eq. JHJ and, especially, of Eq. (jlljl . The 
eigenstates of the Hamiltonian Q can be easily found: 



|0> +a k (E)\k > 



(12) 



i + E*< fcF K(£)l 2 

and the relevant energies are the solutions of the equation 



E 



c/+E 



\Vk\ 2 
E + e k ' 



(13) 



Equation l|13() being superficially similar to Eq. © is to- 
tally different. The former has only real solutions, their 
number beings equal to the number of electron states be- 
low the Fermi level plus one, the latter gives three poles 
(see below), two of which are complex, and the third real 
coincides with the ground state energy E = —Ek- No- 
tice that in the ground state the /-level hole occupation 
number is 



n f 



E 



-E 



n f . 



(14) 



The wave-function right after the measurement can be 
presented as 



m = E 



E 

k<ikp< 



-1/2 



(15) 



Hence at any moment t after the measurement the wave- 
function is 



m = 



E 

E 



E M^l s 



1-1/2 



k<Ck$ 



tp E e' 



-iEt 



(16) 



and the amplitude to find a hole at the level / is 

-l 



a(t) = E 



E 



1+ ]T \a k (E)\ 2 



k<kp 



-iEt 



(17) 



which exactly corresponds to the integral from Eq. 
calculated by the residues method. For t — > oo due to 
dephasing of the contributions from all the states save 
the bound state we obtain 



a(t) 



= n f e lEKt , 



(18) 



coinciding with Eq. flj . 

For a more detailed analysis of Eq. JSJ we additionally 
specify our model, assuming V k = V = const and flat- 
band density of bare itinerant-electron states p = po . In 
this model we get 



E(w) = A In 



D 



(19) 
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where A = |V| 2 po; the imaginary part of the logarithm 
is equal to — n at the real axis between the branch points 
to agree with Eq.Q at the real axis. In this case there 
exist two complex poles. 

Eq. JSJ being valid independently of the strength of the 
hybridization, we however limit our analysis to the Kondo 
regime A -C — e/, where we get 



E K = De~^ A 
E K 



< 1. 



(20) 



One complex pole lo\ — Ef — — ef— in A, with the residue 
approximately equal to 1, is responsible for the tradi- 
tional "Fermi golden rule" decay processes. The second 
pole u>2 = Ek — *ttA with the residue equal to —n f is the 
complex "mirror" of the Kondo pole. As can be shown, 
in the regime we consider, the contribution from the cuts 
is negligible for small and large t. We may hope the at 
the semi-quantitative level it can be neglected at all t, 
though, strictly speaking, the last claim should be sub- 
stantiated by numerical calculations. Thus 



o(i) 



-iE K t--rrAt\ 



,i-7rAt 



(21) 



Due to an entanglement between the localized elec- 
tron and the Fermi sea the local probe of the former 
should create the "decoherence wave" 01 m the conduc- 
tion electron subsystem which will disturb the spin and 
charge distribution in the latter. Let us calculate the 
distribution of conduction electrons around the impurity 
after the measurement. The probability density to find 
the hole at a distance r from the impurity is given by the 
expression 



Pc(r, t) 



\ Tr ikr 1 



, V k e lKr — 

k < kF 2n J -°° u + ek ^ -^- e f 



(22) 



The asymptotic of Ea. (|22|l for large t is obvious: the main 
contribution comes from the real pole u> = —Ek, and we 
obtain 



Pc{r) 



V k e ikr 



k 



Ek - £fc 



(23) 



Notice that the coordinate dependence of the hole den- 
sity of states long after the measurement is the same as 
in the ground state; only the coefficient is smaller by ad- 
ditional factor fif which is just the suppression factor for 
the spectral weight of the Kondo resonance . 

If there exists intermediate asymptotic of exponential 
decay of a(t), it is described by the complex pole Ef. For 
these values of time the probability density is determined 
by the virtually bound state 



Pc{r, t) 



E 



V k e> 



Ikr 



E f - e k 



-2-n-Ai 



(24) 



This asymptotic holds for fif <C e _TrAt <C 1. 

Assuming for simplicity an isotropic hybridization V 
and the dispersion law, we can fulfill angular integration 
in Eas. l|23|l and (|24|l . In the perturbative regime we con- 
sider, Ek is very close to Ep. Due to the logarithmic 
divergence of the integral in Ea. H23|) the main contribu- 
tion comes from the vicinity of the upper limit, and we 
obtain 



Pc{r) 



sm 2 (kpr) _ 2 / Ek — e/ 



(k F r) 
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(25) 



Now we consider another measurement result. Let the 
first measurement has found the localized electron sitting 
at the /-level. The wave-function directly after such mea- 
surement according to the von Neumann "wave-function 
collapse" postulate 0] is again given by Eq. (01 with 
the initial condition 



V>(0)~-/=^a fc /t Cfc!/ |fi> . 

kv 



(26) 



After simple algebra we obtain 



a(t) = — / doj— Y(w), (27) 

2ni J_ 00 E(w) -u-ef 



where 



\v k \ 2 



(uj + e k + is)(E K - £fe) ' 



(28) 



Again closing the contour of the integration and ignoring 
the contribution from the branch cuts we find (in the 
perturbative regime) 



iE K t 




(29) 



^From the orthogonality of eigen- functions corresponding 
to —Ek or ef to the wave- function corresponding to Ek 
we obtain 



Y(-Ek) = Y(e f ) = -1. 
Finally, in the Kondo regime (fif 1) one has 



,iE K t 



Thus for t — > oo 



(30) 
(31) 
(32) 



This result can be understood in a simple way. Since 
in the limit under consideration A <C — e/ the (hole) 
occupation number of the /-level is very small, the mea- 
surement with the result hf = produces negligible dis- 
turbance of the wave-function of the ground state. This 



4 



is in sharp contrast with the measurement which results 
in n/ = 1. After that measurement, as we can see from 
Eq.JTT}, the asymptotic value of nf changes drastically. 

We have analyzed two kinds of processes. The process 
of measurement instantly reduces the ground-state wave 
function, with the probability to find a hole at the /-level 
equal to fif) to the form (jSJ), with the probability to find 
a hole at the /-level equal to 1 (we'll consider here only 
one possible measurement result). After the measure- 
ments starts the process of quantum evolution with the 
characteristic time scale 1/A, which ends up in the state 
with the probability to find a hole at the /-level equal 
to fly. Notice that we obtained seemingly paradoxical 
result. On the physical grounds one should expect the 
disappearance of the measurement effects at large times, 
and hence the return of the hole-occupation number to 
n f . The way to solve the paradox is clear when we look at 
Eq. ifTBl) . Together with the evolution described by this 
equation there exists a third kind of a processes: ther- 
malisation of the hole as the result of the excitation of 
electron-hole pairs. The consideration of such processes 
demands considering instead of the model Nf — > oo the 
model with large but finite Nf and taking into account 
terms of the order of 1/Nf. This will be the subject of 
a separate publication, but already now we can say that 
in such model the time scale of this thermalisation will 
be much larger than the time-scales of the evolution de- 
scribed by Eq. I|21(l and determined by the inverse Kondo 
temperature. Thus even in such model, our results, say 
Eq. Ijl8|) w iH be still valid, only it will present an in- 
termediate asymptotic. Notice that what is said in the 
conclusion is the time-representation of the phenomenon 
of the density of states description of the Kondo reso- 
nance. Whereas in the infinite- Nf limit investigated by 
us the Kondo peak has a zero width, for finite Nf it is a 
Lorentzian with both the distance from the Fermi energy 
and the width of order of the Kondo temperature (and 
for purely spin case Nf — 2 the center of the resonance 
just coincides with the Fermi energy) [3.HH. 

To conclude, we have presented the analytical solution 
of the problem of real time charge dynamics in quantum 
impurity system for the exactly solvable infinite- Nf limit. 
This limit is sufficient to describe the effects of decoher- 



ence on the Kondo resonance. On the other hand, the 
process of recoherence, that is, return to the ground state 
after the measurement, requires the consideration of the 
higher-order processes in 1/Nf. 
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